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Abstract 



In this thesis, we study the critical dynamics near the QCD critical point. 

Near the critical point, the relevant modes for the critical dynamics are 
identified as the hydrodynamic modes. Thus, we first study the linear dynamics 
of them by the relativistic hydrodynamics. 

We show that the thermal diffusion mode is the most relevant mode, whereas 
the sound mode is suppressed around the critical point. We also find that the 
Landau equation, which is believed to be an acausal hydrodynamic equation, 
has no problem to describe slowly varying fluctuations. Moreover, we find that 
the Israel- Stewart equation, which is a causal one, gives the same result as the 
Landau equation gives in the long-wavelength region. 

Next, we study the nonlinear dynamics of the hydrodynamic modes by the 
nonlinear Langevin equation and the dynamic rcnormalization group (RG). In 
the vicinity of the critical point, the usual hydrodynamics breaks down by large 
fluctuations. Thus, we must consider the nonlinear Langevin equation. We 
construct the nonlinear Langevin equation based on the generalized Langevin 
theory. After the construction, we apply the dynamic RG to the Langevin 
equation and derive the RG equation for the transport coefficients. 

We find that the resulting RG equation turns out to be the same as that for 
the liquid-gas critical point except for an insignificant constant. Consequently, 
the bulk viscosity and the thermal conductivity strongly diverge at the critical 
point. Then, a system near the critical point can not be described as a perfect 
fluid by their strong divergences. 

Wc also show that the thermal and viscous diffusion modes cixhibit critical- 
slowing down with the dynamic critical exponents ^thermal ~ 3 and Zviscous ~ 2, 
respectively. In contrast, the sound mode shows critical-speeding up with the 
negative exponent -Zsound ^ —0.8. 
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Chapter 1 

Introduction 



The quantum chromodynamics (QCD) is established as the fundamental theory 
of the strong interaction. Although the fundamental theory is established, we 
can not study a strongly interacting matter based on the first principle at finite 
density. Thus, the QCD phase structure at finite density and temperature is not 
established [U [2] , and determination of the structure is a fundamental problem. 

Figure fTTT] shows a schematic phase diagram of a strongly interacting matter. 
In low density and temperature region, we have the hadronic phase, in which 
quarks and gluons are confined in hadrons. On the other hand, in high density 
and temperature region, the confinement breaks. Then, deconfined quarks and 
gluons become relevant degrees of freedom. The phase boundary line, which 
separate the two phases, is predicted to be a first oder phase transition line by 
various effective models of the QCD [5]. This line exists in the finite density 
region, and then the first principle lattice QCD is not available. In contrast, 
transition along the temperature axis at zero density is predicted to be crossover 
by the finite temperature lattice QCD ^ . Namely, the transition is not associ- 
ated with a thermodynamic singularity. The end point of the first order line is 
considered to be a second order transition point ^ |A^ . This point is called the 
QCD critical point [S]. 

If the critical point exists, the correlation length of an order parameter, ^, 
diverges and thermodynamic quantities have singular behaviors at the point, 
like the specific heat at the liquid-gas critical point. By the singularity, the 
critical point is expected to be useful for experimental probe of the QCD phase 
structure in the relativistic heavy ion collider [B1[71|S]. Thus, the critical point 
attract the interest of many people. 

It is known empirically that, by decreasing colliding energy, the chemical 
potential of the created matter increases. Then, by varying the colliding en- 
ergy, we can experimentally scan a part of the phase diagram. Such experiment 
is called the beam energy scan program and now ongoing j6]. Several experi- 
mental signatures have been suggested based on the critical divergence [71 [TT| . 
For example, the baryon number fluctuation is predicted to be enhanced, as 
{{SN)) near the critical point [l2]. Thus, the baryon number fluctuation is 
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Figure 1.1: A schematic phase diagram of the QCD. 



expected to have the non-monotonic dependence on the coUiding energy, if the 
created matter passes near the critical point. 

The static critical phenomena (namely, time-independent one) have been 
strenuously studied. Consequently, the order parameter, more appropriately, 
the critical mode 9„, is now identified as a linear combination of the chiral 
condensate <j and the baryon number density n |10j . if the critical point ex- 
ists. Moreover, the static universality class is identified as the class of 3d Ising 
model, Z(2)[5]. We note that the chiral symmetry is explicitly broken by the 
finite quark mass, and thus a couples to n. 

In contrast, dynamic critical phenomena (for example, critical-slowing down 
or divergence of transport coefficients) have not been fully studied. The critical- 
slowing down is the phenomenon that the life time of the order parameter di- 
verges at the critical point. The dynamic critical phenomena are of a long-time 
scale. Thus, long-living modes (slow modes) are the relevant modes for the crit- 
ical dvnamicsjl3j. What are the slow modes near the QCD critical ? The slow 
modes are now considered as the hydrodynamic modes coming from the fluctu- 
ations of conserved densities: the baryon number n and the energy-momentum 
7^A»i^ IJOj [J4] _ ^\^(^ slow dynamics of the order parameter, which is the linear com- 
bination of a and n, is governed only by the baron number fluctuation. Thus, 
the chiral condensate would be irrelevant. 

Although the relevant modes are specified, dynamics of them has not been 
studied. Specifically, the coupling between in and 8T^^ is not taken into account 
in the earlier study [14.. Thus, in this thesis, we shall first study the linear 
dynamics of them by the relativistic hydrodynamics. Here, the important point 
is that the long-time behavior of the conserved densities, n and T^^ ^ is basically 
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given by the hydrodynamics. As relativistic hydrodynamic equations, we use 
the Landau equation and the Israel-Stewart equation. The Landau eauation|15) 
is beheved to be an acausal hydrodynamic equation[T6]. However, we show 
that the equation has no problem to describe the hydrodynamic modes. We 
shall also find that the Israel- Stewart equationfl^, which is a causal equation, 
gives the same result as the Landau equation gives on the long-time and long- 
distance scale. Furthermore, we shall show that the actual slow modes is three: 
the thermal, viscous, and sound modes. We also find the relative importance 
of them; the thermal mode is the most relevant, whereas the sound mode is 
suppressed around the critical point. 

Furthermore, some authors suggested a divergence of the bulk viscosity at 
the QCD critical point 18 . But, its validity is controversial |19[ \TT\: for 
example, the limiting operation in the Kubo- formula may not be correct |19] . 
and a study by the relativistic Boltzmann equation |21j shows that the bulk 
viscosity is finite at the critical point. Thus, it is still not obvious whether the 
transport coefficients will diverge or not at the QCD critical point. 

In fact, as is known in condensed matter physics, the critical divergence of 
the transport coefficients is a common phenomenon at a critical point, such as at 
the liquid-gas critical point, and originate from a universal mechanism; nonlinear 
interactions of slow modes cause the divergence \22\ I23j . This implies that fast 
modes, or microscopic processes as described by like the Boltzmann equation, 
would not contribute to the critical divergence of these quantities, if any. The 
dynamic renormalization group (RG) theory |241I13| is a standard technique for 
the critical dynamics, which systematically incorporate the macroscopic nonlin- 
ear interaction causing the divergence of transport coefficients. In this theory. 
We must construct a nonlinear Langevin equation to describe the nonlinear 
interaction of the slow modes. The construction goes as follows. First, We iden- 
tify the slow variables, which label a state on the long-time and long-distance 
scale. Next, the thermodynamic potential for the slow variables is constructed 
to determine the static property of the system. Finally, the streaming terms, 
which cause the dynamic-nonlinear interactions, and the kinetic coefficients are 
determined for respectively describing time reversible and irreversible changes 
of the slow variables. We note that the streaming term is absent in the simple 
Brownian motion. 

The general theory of the critical dynamics as described above tells us that 
an essential ingredient is to properly construct the nonlinear Langevin equation 
for the critical dynamics. As far as we know, this is the first attempt for the 
QCD critical point. Our construction of the Langevin equation is based on the 
generalized Langevin theory, so-called the Mori theory [351 US] ; and the relativis- 
tic hydrodynamics, because the slow modes are identified as the hydrodynamics 
modes |10|, I14|, I26j ; we construct the streaming terms from continuity equations 
and the potential condition, which is a general condition for streaming terms 
[24j [13]. Also, we use the thermodynamic potential for the 3d Ising system 
as that for the QCD critical point because the static universality class is the 
same as 3d Ising class [lOl [Mj [27] . Finally, we determine the kinetic coefficients 
from a relativistic hydrodynamic equation, here the Landau equation [15) used. 
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Consequently, we shall show that the Langevin equation differs from it for the 
liquid-gas critical point by relativistic effects, although the dynamic universality 
class of the QCD critical point is conjectured as of the liquid-gas critical point 

After such construction, we apply the dynamic RG to the Langevin equation 
and derive the RG equations for the transport coefficients. Consequently, to our 
surprise, these RG equations turn out to be the same as for the liquid-gas critical 
point except for a irrelevant constant, although the Langevin equations are 
different. Therefore, the bulk viscosity and the thermal conductivity strongly 
diverge and can be more important than the shear viscosity near the QCD 
critical point. We shall also show that the thermal and viscous diffusion modes 
exhibits critical slowing down, whereas the sound mode critical speeding up. 

This thesis is organized as follows. 

In Chaplll we review the general theory of critical dynamics. This thesis is 
based on the general theory. Specifically, we first give the projection operator 
method. By this method, we can systematically decompose any dynamic vari- 
ables into a slowly varying motion and a rapid one. Namely, we can extract the 
relevant motion for the critical dynamics. 

We also give the statistical basis of coarse-grained equations of motion, 
namely, the Langevin equation and the hydrodynamic equation. We note that a 
microscopic theory is difficult to describe the long-time effect although the dy- 
namic critical phenomena are of the long-time scale. Thus, the coarse-grained 
equation is relevant. 

After the projection operator method, we briefly explain the typical dynamic 
critical phenomena: the critical-slowing down and the critical divergence of 
transport coefficients. We also give the concept of the dynamic RG, and the 
critical and hydrodynamic regimes. Near the critical point, the macroscopic 
scale is divided into the two regimes. 

Moreover, we also give the earlier studies on the slow modes near the QCD 
critical point. 

In ChaplHl we study the linear dynamics of the slow modes by the relativistic 
hydrodynamics . 

In ChapUl we study the nonlinear dynamics of them by the dynamic RG. 
In the final chapter, we give summary and concluding remarks. 
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Chapter 2 

Theory of critical dynamics 



Here, we provide the general theory of critical dynamics, which is developed in 
condensed matter physics. 

We also give earlier studies [TOl [H] on the QCD critical point, which shows 
the relevant variables are the hydrodynamic ones. 

2.1 Projection operator method 

The Mori's projection operator method gives the microscopic basis of a Langevin 
equation |251I23| . By this method, we can formally extract dynamics in the long- 
time scale, which is relevant for critical phenomena. The Langevin equation has 
been widely used to study critical dynamics. 

2.1.1 Linear Langevin equations 

Here, we derive a linear Langevin equation from a microscopic equation by the 
Mori theory [m [53. 

Now, let us consider a classical many body system, for simplicity. A general- 
ization to a quantum system is straightforward. For the classical system, a time 
evolution of an arbitrary dynamic variable is given by the Liouville equation : 

^^B{t)^{B{t),H}pB., (2.1) 

where B{t) is a dynamic variable in time t, H a microscopic Hamiltonian and 
{, }pB the Poison bracket. Introducing the Liouville operator, iC, as 

iCB{t) = {Bit),H}pB- (2.2) 

we can formally solve Ea. (|2.ip as 

B{t) ^ e'^^B{t ^0). (2.3) 

^If we replace the Liouville equation with the Hisenberg equation, the following arguments 
are valid for a quantum system. 
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In the following, we shall decompose the time evolution Eq.(| 

First, let us represent a set of slowly varying variables (slow variables) as 

{Aj{t)}. The slow variables label a macroscopic state and describe a time 

evolution on a macroscopic scale. 

Next, we define the linear projection operator V as 

rB{t)=Y,{B{t)A,)xJ,'A,. (2.4) 

jk 

Here, (...) is the equilibrium-statistical average, Aj an initial value of the slow 
variable, namely Aj = Aj(t — 0), and xjk inverse of the correlation Xjk = 
(AjAk). In the following, we denote the initial values of the slow variables 
without the argument t. 

The operator V extracts a slowly varying motion from arbitrary dynamical 
variables. Also, we define the orthogonal operator, as Q = 1 — "P. 

Now, we use the operator identity, which is valid for arbitrary iC and V |24j . 

-^e'^* = e*^*Pz£ + / dte'^'^*~''^PiCe^'^''Qi£ + e^'^'QiC. (2.5) 
ot Jo 

Multiplying Ea. (|2.5l) by the initial values of the slow variables Aj, we obtain 
the linear Langevin equation for Aj{t) — exjp[iCt\Aj: 



(t-t) + f,{t), (2.6) 



k k 

without any approximations. Here, we introduced 

^Y.^A,Ai)xik\ (2.7) 
fe 



fj{t) =exp(Qi/:t)Qi,, (2.8) 

^3k{t) =Y.^fAi)m)xTk'^ (2.9) 
I 

with Aj = iCAj. The equation (|2.6p has the following properties. 

1. Eq. (|2.6p is the exact relation. Here, we only used the operator identity. 

2. The first term in the right-hand side is a time-reversible change. 

3. The second term is a time-irreversible change. Also, this term depends on 
a past time value, Ak{t — i ). Tjk[k) is called a memory function. 

4. The last term is a rapid motion and usually treated as a random noise. 

Now, we give a transport coefficient in this scheme. If the time-scale of 
the slow variables Aj{t) and that of the noise fj{t) are well separated, we can 
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assume that Ak {t—t ) does not change among the correlation time of the memory 
function, Eq. (j2.9p . Namely, we approximate the time- irreversible term as 

dtTjk{t)Ak{t-t) ^ (^j^ dtT,k{t)^Ak{t). (2.10) 

This approximation is called the Markov approximation. We now introduce the 
linear transport coefhcient as 

" dtV,i{t)^xik, (2.11) 

{fAt)hm- (2.12) 

We see that the transport coefficient is given as the time correlation of the 
noises. Finally, we obtain 

^^A,{t) ^Y.'^l,Mt) -Y.h^^kMt) + f,{t). (2.13) 

k hi 

We see that Ea. (|2.13p loses the memory effect. 

An important point is that the noise, fj{t), implicitly includes nonlinear 
terms of the slow variables. Namely, fj{t) is not orthogonal to the nonlinear 
terms: 

ifjm,) =0, (2.14) 
{f,it)AkAi) ^0. (2.15) 

This originates from that T' is a linear projection operator. Therefore, if we can 
not neglect the nonlinearity, fj{t) can not be treated as a noise. 

Moreover, from Eq. (j2.12p . we see that the nonlinear terms contribute to the 
linear transport coefficients. This contribution causes the critical divergence of 
the transport coefficients near a critical point. Near the critical point, fluctua- 
tions become large. Therefore, we can not neglect the nonlinear fluctuations and 
must consider a nonlinear Langevin equations as a basic equation for critical 
dynamics. 

2.1.2 Nonlinear Langevin equations 

Here, we derive the nonlinear Langevin equation by the nonlinear projection 
operator [Hdl]. 

To define the nonlinear projection operator, we first introduce the following 
delta functional as 

g{A,a)^l[SiA,~a,), (2.16) 

3 

where aj are some initial values. The equilibrium-statistical average of this gives 
the equilibrium-distribution function: 

P,^{a)^{g{A,a)). (2.17) 
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With this delta functional, we can define the nonlinear projection operator 
acting on any dynamic variables B as 

rniB = {Bg{A,a))/P,^{a). (2.18) 

The physical meaning is simple. We fix the slow variables Aj at some values aj 
and average out the other degree of freedom. In other words, we eliminate the 
fast variables and extract a slowly varying part that is determined by only the 
slow variables. 

The important point is that the nonlinear projection on Aj is identical to 
the linear projection on g{A, a). The linear projection on g{A, a), which is Vg, 
is given by 



VgB = J dada {BgiA,a)){g{A,a)giA,a))-'giA,a), (2.19) 

=.{Bg{A,a))/P,^{a), (2.20) 
= VniB. (2.21) 

Here, we used the relations 

{g{A,a)g{A,a)) = ,5(a - a')Peq(a), (2.22) 
{giA,a)giA,a))-' = J(a - a )/^eq(a). (2.23) 

Namely, the linear projection Vg is equivalent to the nonlinear projection "Pni- 
Hence, we can derive the nonlinear Langevin equation about Aj from the linear 
Langevin equation about g{A, a). However, we leave the derivation to Appendix 

El 

The resulting nonlinear Langevin equation with the Markov approximation 

is mm] 

k 

with /5 being the inverse temperature. Here, we introduced 

v,{a) =(i,;a), (2.25) 
0j{t) ^exp[Qgia]QgAj{0), (2.26) 

Uk{a) = / dt{e,{t)9k{Q)\a), (2.27) 



where (...; a) — {...g{A^ a))/Pcq{a) is the conditional average in which Aj is fixed 
at aj. Also, we defined the thermodynamic potential (or the effective potential) 
H{A) as 

1 

Z 

where Z is a normalization constant. 

Now, we give physical meanings of Eq. (|2.24p . 



Peq(A) = -cxp[-/3iJ(A)], (2.28) 
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1. The first and second terms are the slow motions and nonlinear in Aj. 

2. The first term, which is called the streaming term, gives a time-reversible 
change. 

3. The second term gives a time-irreversible change. Lik{a) is called the bare 
transport coefficient. 

4. The last term is a fast motion and treated as a stochastic variable obeying 
the fluctuation-dissipation relation 

{ej{t)ek{t);a) ^2L^k{a)5{t~t)- (2.29) 

In contrast to the linear case, Oj{t) does not include the slow variables. 
Namely, the nonlinear terms of Aj are explicitly extracted in the first and 
second terms. 

Even for the QCD critical point, we may use the generalized Langevin equation, 
because only the time-scale separation is assumed in the Mori theory. Further- 
more, we note that, by the time-scale separation, transport coefficients arises 



2.1.3 On slow variables 

In the Mori theory, a choice of the slow variables plays a crucial role. How 
we choose the slow variables? On a macroscopic scale, the slow variables are 
given as conserved densities, Nambu-Goldstone modes (NG modes), and order 
parameters [13] . We now explain why they are slow. 

First, let us consider the conserved density. The important point is that any 
conserved densities generally obey a continuity equation; 

^ = -V.i(r,.), (2.30) 

where n is a conserved charge density and j is its current density. Performing 
Fourier transformation about r, we have 

^^ = -^k■jik,t), (2.31) 

where k is the wavenumber. We see that the time-change rate is proportional 
to the wavenumber. Thus, the conserved density is slow in the low-wavenumber 
region, namely, the macroscopic scale. 

Next, we consider the NG modes. The key is that the NG modes gener- 
ally have gapless-dispersion relations. Namely, their dispersion relations are 
proportional to the wavenumber or the square of that: 

uj{k) (xk or k^, (2.32) 

where a; is a frequency. Again, in the low-wavenumber region, we have the slow 
motion for the NG mode. 
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Finally, let us consider the order parameters. In general, near a critical 
point, the relaxation of the order parameter is anomalously slow. Such the 
dynamic critical phenomenon is called the critical-slowing down. Thus, the 
order parameter is slow near the critical point. We shall explain about the 
critical-slowing down in the next section. 

Now, the important point on the slow variables is that the low-wavenumber 
components are slow; but the high-wavenumber components are fast. Thus, we 
must restrict the wavenumber by the ultraviolet cutoff, A. Then, the Langevin 
equation has the ultraviolet cutoff. 



2.2 Critical slowing down 

Here, let us illustrate the critical-slowing down by a linear Langevin equation. 
For example, we consider the case that a single-order parameter, (f), is only the 
slow variable. In this case, we have the following linear Langevin equation: 

^ = ^Ht) + m. (2.33) 

Here, a reversible term is absent. The reason is the following. From Ea. (|2.7p . 
we have the relation ^ISj , 

in,, = ({A„Aj}pb) or {[A,,A,]/{ih)). (2.34) 

Here, [..., ...]/{ifi) denotes commutation relation for a quantum system. There- 
fore, if we have a single slow variable, the reversible term is generally absent. 

Now, the important point is that the susceptibility of the order parameter, 
Xipipi generally diverges at the critical point. Thus, the relaxation of the order 
parameter exhibits slowing down near the critical poinlH. 



2.3 Hydrodynamic and critical regimes 

Here, we give a valid region of the linear and nonlinear Langevin equations in 
terms of the wavenumber and the cutoff [24l [28] . 

Now, let us consider a non-equilibrium state near an equilibrium state. In 
other words, the state fluctuating from the equilibrium state is considered. We 
note that the critical point is defined on the equilibrium phase diagram. For such 
state, the dynamic variables turn out to be fluctuations from the equilibrium 
state. Then, if we consider a state far from the critical point, the fluctuations 
are small and the linear Langevin equation suffices. 

In terms of the wavenumber, the valid region is 

0<k^a-\ (2.35) 

^ As we shall see in the sectioi l274l the transport coefficient, Lif,^, also diverges. However, 
the divergence of the transport coefficient is typically weaker than that of the susceptibility. 
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where a is a microscopic characteristic length scale. Thus, the ultraviolet cutoff 
of the linear Langevin equation, Al, is chosen as 

Al < a-\ (2.36) 

in the normal region. 

In contrast to the normal region, near the critical point, the fluctuations 
become large. Thus, we can not neglect the nonlinear fluctuations and must ba- 
sically use the nonlinear Langevin equation. However, if we restrict our interest 
to a much larger scale than the correlation length of the order parameter, ^, we 
can again use the linear Langevin equation even in the critical region. On such 
scale, the information of the critical point is obscure and included in parameters, 
like the transport coefficients. Namely, the cutoff for the linear theory must be 
chosen as 

AL«r\ (2.37) 

in the critical region. 

In contrast, the nonlinear Langevin equation can describe the nonlinear fluc- 
tuations and thus is valid even on a smaller scale than ^. Then , we can choose 
the cutoff for the nonlinear Langevin equation as 

Anl < (2.38) 

even in the critical region. 

The important point is that, near the critical point, the wavenumber regime 
is divided into the two regimes: hydrodynamic and critical regimes. 

The hydrodynamic regime is 

0<fc<r^. (2.39) 

In this regime, the linear theory is still valid E|. On the other hand, the critical 
regime is 

< < (2.40) 

In this regime, only the nonlinear theory is valid and fully reflects the informa- 
tion on the critical point. 

We shall discuss the two regimes in terms of the dynamic RG in Sec l2.5l 

2.4 Critical divergences of transport coefficients 

The critical divergence of transport coefficients (or diffusion constants) is a 
common phenomenon, for instance, at the liquid-gas critical point, ferromag- 
netic transitions and so on [24l [13]. The important point is that the critical 
divergence originates from a universal mechanism; macroscopic nonlinear fluc- 
tuations, namely, the nonlinear terms of slow variables, cause the divergence, 
which is implied in Sec I^XT] [H US]. 

^The linear Langevin equation is sometimes called hydrodynamic theory. So, this regime 
is called the hydrodynamic regime. 
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Now, we illustrate how the macroscopic nonlinear fluctuations cause the 
critical divergence. For an example, let us consider the thermal conductivity 
near the liquid-gas critical point [29] . The thermal conductivity is given by the 
Kubo formula as follows, 

X = T-^JdrJ^ dt{q{r,t)q{0,0)), (2.41) 

where q{r,t) and T are the heat current and temperature, respectively. The 
heat current q{r, t) consists of two parts: one is due to a microscopic process as 
described by a microscopic theory, like the Boltzmann equation, and the other 
is by the nonlinear fluctuations of macroscopic variables |25j : 

Q ^micro ^macroi (2.42) 

where (/micro and gmacro respectively denote the microscopic and macroscopic 
currents. The macroscopic process causing the heat current is identified as the 
entropy density convected by fluid velocity fluctuation. Thus, we have 

ymacro 

6sSv, (2.43) 

where Ss and Sv respectively denote the fluctuations of the entropy density 
and the fluid velocity. The macroscopic current Eq. (I2.43P is of the second 
order in fluctuations and hence negligible far from the critical point. However, 
it becomes the main part near the critical point, because the fluctuations are 
enhanced there. We see that Eq. p.41|) now has the following form 

A - Amicro + y dt{Ss{r,t)Sv{r,t)Ss{0,0)Sv{0,0)), (2.44) 

where Amicro is the thermal conductivity coming from qmicro- Recalling that 
the entropy density fluctuation is the order parameter for the liquid-gas critical 
point, we see that the second term of Eq. (|2.44p diverges at the critical point. 
This is the mechanism causing the critical divergence of transport coefficients. 

Let us call the transport coefficients, such as Amicro 7 coming from microscopic 
processes the bare transport coefficients, and those including the contributions 
from the nonlinear macroscopic fluctuations the renormalized ones. 

Then, we need not to study the critical divergence of transport coefficients 
by a microscopic theory because the divergence originates from only the macro- 
scopic processes. The dynamic RG [Ml lEl |3TJ [32l is the standard theory 
treating such nonlinear macroscopic fluctuations. In this theory, we must con- 
struct a nonlinear Langevin equation as a basic equation for the critical dynam- 
ics. 

Furthermore, we note that the earlier studies [13 [21], on the transport co- 
efficients near the QCD critical point, treat the bare part. Thus, those do not 
take into account the contribution from nonlinear macroscopic fluctuations. In 
Chaplll we shall study the renormalized part by the dynamic RG. 
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2.5 Dynamic RG 

Here, we give a conceptual aspect of the dynamic RG. A technical aspect is well 
given in the textbook [13] . 

The general dynamic RG transformation usually consists of two procedures, 
i.e., coarse graining and rescaling as in the static RG transformation [30l [13] . 
However, as is shown in [2H[17J[11], we can omit the rescaling, if we are interested 
in only the critical exponents of transport coefficients. 

The nonlinear Langevin has a ultraviolet cutoff Ag, which should satisfy the 
following inequality 

r^<Ao<a^^ (2.45) 

Namely, at the starting point of the dynamic RG, we are in the critical regime. 

Then, the Langevin equation is coarse grained by averaging over the high- 
wavenumber components of the slow variables Aj (t) in the infinitesimal wavenum- 
ber shell, 

A - (5A < fc < A, (2.46) 

for Eq. (|2.24p . Here, A starts from the initial value Aq and is lowered up to 
A <C Namely, at the final point, we are in the hydrodynamic regime. The 
eliminated components turn out to be included in parameters of the Langevin 
equation. Thus, the parameters, like the transport coefficients, are renormalized. 

In other words, we first construct the nonlinear Langevin equation to de- 
scribe nonlinear effects in the critical regime, and the nonlinear effects are in- 
cluded in the liner transport coefficients in the hydrodynamic regime by the RG 
transformation. 

2.5.1 Contrast with the static RG 

Here, we first stress that the concept of the dynamic universality class is not so 
universal contrary to its name. Then, the class of the QCD criticalpoint may 
not be the same as of the liquid-gas critical point or the model H □, although 
it is conjectured by [HI [20]. To see this, let us contrast the difference between 
the static RG with the dynamic one. 

An important point is that the respective infrared effective theories are dif- 
ferent; in the static case, the infrared effective theory is the thermodynamic 
potential (or so-called Landau free energy). Then, its function form about or- 
der parameters are determined only by the space dimension and the symmetry 
among the order parameters but not by microscopic details . Thus, the concept 
of the universality class makes sense for the static case. In contrast, for the 
dynamic RG, the infrared effective theory is the nonlinear Langevin equation. 

Here, the important difference arises; the relevant variables for the Langevin 
equation is not only the order parameters but also conserved densities and NG 

*The model H 1221 1281 is the minimal-dynamic model for a critical point that its relevant 
modes are given as the nonrelativistic-hydrodynamic modes. The liquid-gas critical point 
belongs to the dynamic universality class of the model H 
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modes, and its nonlinear couplings can not be determined by only the symme- 
try in general. Consequently, the dynamic universality class is not so universal 
compared to the static one. Specifically, the nonlinear couplings, namely, the 
streaming terms vj (A) , are generally given by the Poisson brackets (commuta- 
tion relations) among the slow variables in the classical (quantum) system; 

k 

where 

QjkiA) = {{Aj,Ak}pB;A) or {[A,, Ak]/{ihy, A). (2.48) 

The above expression, Ea. (|2.47l) . is derived from Eg. (12. 251) . see [13] for the 
derivation. The important point is that the Poisson-bracket relations depend 
on the microscopic expressions of the variables. This fact leads to an important 
consequence that the dynamic universality class of the QCD critical point may 
not be the same as of the liquid-gas critical point or the model H. Actually, in the 
model H, the Poisson-bracket relations are calculated with the non-relativistic 
relations [22l[T3]. 

2.6 Slow variables near the QCD critical point 

Here, we give earlier studies on the slow variables near the QCD critical point. 
The slow variables have been identified as the fluctuations of the conserved 
densities: the baryon number, n, and the energy and momentum, T'"^. 

First, let us consider slow variables for chiral limit although we are interested 
in the finite quark mass case. For the chiral limit, the slow variables are a, tt, 
and the conserved densities. Here, a is the order parameter about the chiral 
phase transition and tt is the NG mode for the spontaneous chiral-symmetry 
breaking. Thus, the above quantities are slow. The critical dynamics in this 
case is studied in [39l |40] . 

In contrast to the chiral limit, for finite quark mass, the chiral symmetry is 
explicitly broken. Thus, the pion has mass and is fast. Moreover, by the explicit 
symmetry breaking, a mixes to the conserved densities. H. Fujii and M. Ohtani 
showed that, by this mixing, a also becomes massive, and the fiat direction of 
the thermodynamic potential is a linear combination of a and n [lOj . 

Figure \2A\ which is based on NJL model and adapted from [10], shows the 
thermodynamic potential in (cr, n) and (ct, s) planes, where, s is the entropy 
density. Here, (a) and (c) respectively denote the potentials for the chiral limits 
and for the finite quark mass. Although (b) denotes for a tricritical point, we 
do not treat the point, see [10] for the detail. We see that, for (a), the potential 
flats along the a direction whereas, for (c), the potential flats along the linear 
combination of a and, n or s. We note that, in this study, the quark mass is 
only a few MeV. 

Moreover, D.T. Son and M.A. Stephanov showed, by a linear Langevin equa- 
tion, that the long-time behavior of the linear-combination mode is determined 



SH 
5Ak 



SA, 



(2.47) 
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Figure 2.1: The thermodynamic potential in {a, p) and (a, s) planes. Here, J7 is 
the thermodynamic potential, p the baryon-number density, and s the entropy 
density. For (a), the potential flats along the a direction whereas, for (c), the 
potential flats along the linear combination of a and, p or s. This figure is 
adapted from [TU] . 



only by the conserved densities [H]. Thus, a mode just traces the conserved 
densities and is unimportant. 

Now, we briefly give the study, [M]. First, to connect our Langevin equation, 
Eq. (|2l^ . to that in [M], we slightly rewrite Eq. ([2J3l) as 



^^A,{t) = Y,^n,Mt)-Y,LJ-^^^^ + Mt). (2.49) 

fe k 



Here, we introduced the thermodynamic potential as the Gaussian form: 



(2.50) 



If we substitute this potential in Eq. (|2.49p . we get back Eq. (|2.13p . In the fol- 
lowing, we consider Eas. (|2.49p and (|2.50l) . 
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For simplicity, we now neglect the energy and momentum densities and 
consider the mixing only between the fluctuations of a and n. Namely, our 
slow variables are 

{Aj} = {5(j,6n} (2.51) 

In this case, the potential i^ 



PH{5u, 5n) = J dr 



(2.52) 



where A, B and C is related to inverses of the susceptibility. The important 
point is that the second term explicitly break the chiral symmetry: a — > —a. 
By this term, we take into account the quark mass. 

Then, we have the linear Langevin equation for a and n in the Fourier space: 

at on OCT 

Here, time- reversible terms are absent by Eg. (12.341) and the time-reversal sym- 
metry. In general, those terms vanish if we have only variables whose time- 
reversal properties are even. 

Now, let us consider the transport coefficients. First, from the Onsager's 
reciprocal relational], we have 

L„,{k) = Lnaik). (2.55) 

Furthermore, we now expand the wavenumber dependence 

L„,{k) ^T + 0{k^), (2.56) 
L,„(fe) ^Xk^ + 0{k^), (2.57) 
Lnn{k) ^Xk^ + 0{k^), (2.58) 

because we are interested in the low- wavenumber region. Here, F, A and A are 
wavenumber- independent constants. The important point is that the expansions 
of L„cr and L„„ start at the order k^. The reason is that the baryon number 
density n is conserved, and thus its time-change rate must vanish in the limit 
fc ^ 0, see Eqs.(I23I|) and (P3i|) . 

Then, we finally arrive the linear Langevin equation at the leading order in 

k: 

^^TASa-TBSn + U, (2.59) 
= -{\A + XB)P6a - {XB + XC)k^Sn + /„. (2.60) 



dt 



'In 1141 . derivative terms are included, but these are needless. 
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From the above equations, we can obtain its eigen modes and dispersion rela- 
tions. 

Now, we leave the derivation to [14] and give only the results. The eigen 
modes are the following two modes: 

- BSa + ASn and Sa. (2.61) 

We see that these are the linear combination of a and n, and only a. These 
dispersion relations are 

^a+n{k) ^-tDk\ (2.62) 

LOa{k) ^ -iTA + 0{k'^), (2.63) 

where we introduced the diffusion constant as 

D = XC~X—. (2.64) 

We see that the linear-combination mode is slow, whereas the sigma mode has 
the gap in its dispersion, and thus it is fast. Namely, the slow dynamics of 
the linear-combination mode is determined only by the density fluctuation 6n, 
whereas Sa just traces Sn. Then, we have the only one slow mode. 

Now, we also write the spectrum (time correlation in Fourier space) of them: 

for a later comparison. 

In this study, the coupling with the energy and momentum is not taken into 
account. In the next chapter, we shall consider the couplings. As a result, we 
shall find that actual slow modes is three: thermal, viscous and sound modes. 
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Chapter 3 



Linear dynamics of the 
hydrodynamic modes by 
relativistic hydrodynamics 



In this chapter, we study the slow dynamics near the QCD critical point in the 
hydrodynamic regime. 

The slow variables for the QCD critical point is identified as the fluctuations 
of the conserved densities, as shown in the previous section. Thus, we have the 
slow variables: 



In a straightforward way, we must construct the linear Langevin equation 
for them from Eqs. (|2.13p and ()2.34p . However, we here develop the linear equa- 
tion by linearizing relativistic hydrodynamic equation. We note that the slow 
dynamics of the conserved densities is basically given by the hydrodynamics. 
Then, the resulting linear equation is equivalent to that by the straightforward 
derivation. 

As relativistic hydrodynamic equations, we use the Landau equation and 
the Israel-Stewart equation. The Landau equation is believed to be an acausal. 
Namely, a propagation speed of an information in the equation is considered 
to be faster than the light speed. In contrast, the Israel-Stewart equation has 
relaxation times, and the causality problem is formally resolved. 

However, we shall show that the Landau equation has no problem to describe 
slowly varying fluctuations. Furthermore, we shall find that the Israel-Stewart 
equation gives the same result the Landau equation gives in the long- wavelength 
region. 

The relativistic hydrodynamic equation is given by the following conservation 
laws 



{Aj} = {Sn,de^{ST' 



),SJ' ^ {6T°'')}. 



(3.1) 



d^N^ - 0, 



(3.2) 
(3.3) 
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where N'^ and T^'' are the baryon-number current and the energy-momentum 
tensor, respectively. Those are given as 

N^' = Tiu^ + i/^, (3.4) 
Tt"" = hu^'u" ~ Pg^'" + T^"", (3.5) 

where h = e + P is the enthalpy density with e and P being the energy density 
and the pressure. Also, = (7, 71;) are the fluid four velocity, with 7 being the 
Lorentz factor, and the dissipative terms, i'^ and r^''. The dissipative terms 
differ among considered equations. 

3.1 For Landau equation 

For the Landau equation, the dissipative terms are 



^X{ — ]d'im, (3.6) 

CA^'id^-u), (3.7) 



where A, ry and C are the thermal conductivity, the share and bulk viscosities, 
respectively. A^'^ = g^^ — u^u'^ is the projection onto the space-like vector and 
d'^ = A^'^di, is the space-like derivative. 

Now, we linearize the Landau equation about fluctuations from the equi- 
librium values. Let us write n{x) = ric + Sn{x), e{x) = ec + Se{x), P{x) = 
Pc + SP{x), {I3fi){x) = {I3fi)c + S{l3fi){x), and u'^(a;) = < + (5u^(x). Here, 
the symbols with a preflx S denote the fluctuations. The equilibrium values 
are denoted by a suffix c. Hereafter, variables with the suffix and the preflx 
respectively denote the equilibrium values and fluctuations. 

For simplicity, let us choose the rest frame as the reference frame: = 
(1,0). Then, by the relation u'^Su^i = 0, we have the fluid-velocity fluctuation 
as 

Su^' = {0,Sv). (3.8) 

We also note that the fluid-velocity fluctuation is related to the momentum 
density as 

SJ = hcdv (3.9) 
Then, Landau equation, Eas. (|3.2l) - (|3.7p . are linearized as 

dSn ^ c . \ ( ^^cTc^ ^ 

^= -KV-Sv, (3.11) 
ot 

^ = -V(<5P) + Uo + l^Vo) V(V • 6v) + i^aVHv. (3.12) 



n,V-(5t; + Ao( ^ ) V^<5(/3m), (3.10) 
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This is the hnear equation of motion for our slow variables 0. This equation 
includes the couphngs among the baryon number and the energy-momentum. 
In the following, we study effects of the couphngs on the density fluctuation 
Sn. Namely, we shall derive a spectral function of the density fluctuation from 
Eqs. (j3A0| - ([3T^ . 

Now, we have five equations for seven unknown quantities, Sn, 5e, SJ, 5P, 
and S{/3fi). To solve these equations, let us expand the thermodynamic quanti- 
ties, Se, SP, and 5{I3^), with the density and temperature fluctuations: 

de^ 5n+(^ dT, (3.13) 



The merit of the set {6n, ST) is that their equal-time correlation is orthogonal 

{6nik,t)dT{-k,t)) =0, (3.16) 

for a grand canonical ensemble [231 US] ■ 

In terms of {dn,6T,Sv), the equations (I3.10p - p.l2p take the form 

— —A— V dn + ?ioV • 

ot hcCp J 

+x!^(l-±£l^)wHT = 0, (3.17) 
hc\ Cp J 



dt ' V 3 

hcC^c 
noCp 



+^^^ Vfa + ^^^'^'-"^ ^ 0, (3.18) 



hcC^Cvap d 



noCp dt no7 



naCv d ^^f clcyap 1 



^ . rj. >■ 5T^Q, (3.19) 

Jo ot \ Cp io/ 

where c„ — To{ds/dT)n and Cp — To{ds/dT)p are the specific heats at constant 

1/2 

volume and pressure, respectively, Cg = {dP/de)a the sound velocity, ap — 
— {\/nc){dn/dT)p the thermal expansivity at constant pressure. Here, we used 
some thermodynamic identities, see ^26j for detail. 

Now, let us perform Fourier-Laplace transformation, like 

/+00 I'ca 
dr / dte-^*-'^-^(5n(r,t). 
-00 Jq 



^Wo can show that Eas. l|3.10| l- I|3.12| l are identical to the linear Langevin equations derived 
from Eqs. 1 12. 131 1 and 112.3411 . straightforwardly [38) . 
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Then, we find 



flcCp 



5n{k^ z) + iuck • Sv{k^ z) 



+k^X^ ( ^Sn^fll _ 1 ) - 5n{k, t = 0), (3.20) 



1 



zhc + k^T] ] Sv{k, z) oV] ■ 5v{k, z)) + ik " 6n{k, z 



+ik^^^^^^^5T{k, z) = hcSv{k, t = 0), (3.21) 



cp 



71 (^Cp 



n^Cp 



5n{k, z) + 



clcyap 1 



ST{k,z) 



^^^^^^5n{k,t = 0)+'^ST{k,0). (3.22) 



Here, we note that the mitial values, like dn{k,t — 0), arise from the time 
derivative terms, because we performed the Laplace transformation about time. 

It is convenient to divide the velocity into longitudinal and transverse com- 
ponents 

Sv\i{k,z) EE k ■ Sv{k,z) (3.23) 
6v±{k,z) = 6v{k,z) ~ kSv\i{k,z) (3.24) 

The transverse component of Eas. p.20p - (l3.22p reads 

{zhc + k^rf)5v^{k, z) = hc6v±{k, t = 0). (3.25) 

Now, let us first study the transverse component. The solution is given by 

fa±(fc,0) 

'""^''^^ = . + (WMfc^ - 

Performing the inverse Laplace transformation 

1 



(3.26) 



Sv±{k,t) = 



2TTi 



(5— zoo 



dze Sv±{k, z). 



we have 



SvA^{k,t) = e~(''/'''=)'='*(5t;_L(fc,0). 



(3.27) 



(3.28) 



We see that the transverse component of the momentum diffuses, without prop- 
agation. This modes is called the viscous diffusion mode. 

Furthermore, let us derive spectral function of Svj_. Performing Fourier 
transformation about time t, we obtain 



Sv±{k, uj) 



w2 + {ri/hcYk^ 
28 



5v^{k,Q). 



(3.29) 



3.1. FOR LANDAU EQUATION 



Multiplying Ea. p.29p by the initial value dv±{k,0) and taking the statistical 
average, we have the spectral function 



{{Sv^{k,0)y). 



(3.30) 



Here, we note that {{Sv±{k, 0))^) does not have information on the time evolu- 
tion. 

Now, Let us return to the longitudinal component. The longitudinal com- 
ponent of Eqs. p.20p - p.22|) can be written as the following matrix form 



(Sn{k, z) ' 
dv\\ik,z) 
ST{k,z)^ 



6n{k,0) \ 
(5u||(fc,0) 
-Sn{k,0) + n^JT{k,0)j 



(3.31) 



where the matrix A is 



ik- 



ikric 
z + vik"^ 



TlQCp 



"ha 
ih 

L2Cp 



(3.32) 

Here, we introduced the longitudinal kinetic- viscosity j//, and the thermal diffu- 
sion constant D^, 



UcC. 



(3.33) 
(3.34) 



Multiplying the inverse A ^ from the left in Eq. (|3.3ip . we obtain the Fourier- 
Laplace coefficient of the density fluctuation 



Sn{k, z) 



{A-') 



13 



5n{k,0) + {A-')i25viiik,0) 



+ ^(A"i)i3<5r(fe,0). (3.35) 



Here, an important point is that 6n is orthogonal to ST and Sv: 

{6n{k,0)ST{k,0)) =0, (3.36) 
{Sn{k,0)6vii{k,0)) =0. (3.37) 

The second equation comes from the time-reversal invariance of the equilibrium 
state. 
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Thus, by the similar procedure as in the transverse component, we obtain 
the spectral function of the density fluctuation 

Snn{k,uj) = {Sn{k,uj)5n{k,t = 0)} 



{{Sn{k,t = 0)f) 



j2 + D^k^ 



where we have introduced the sound diffusion constant, Z^s, as 

r.2 



1 



1^1 



2D,ap 



(3.38) 



(3.39) 



The detailed derivation is given in Appcndix[Bj We see that the spectral function 
has three peaks at frequencies uj = and cj = ±c<;fc: The peak at a; = 
corresponds to thermally induced density fluctuations. This mode is called 
the thermal diffusion mode. The two side peaks at w = ±Csk correspond to 
mechanically induced density fluctuation, i.e. sound waves. This mode is called 
the sound mode. Roughly speaking, if we expand the density, Sn, with the 
entropy density, Ss, and the pressure, SP, 



dn\ 



Sn = [ — ] 5s + [ — ] 5P, 



ds 



dn 



dP 



(3.40) 



we see that the first term corresponds to the thermal mode, while the second 
term the sound mode. 

Now, let us compare this result with that in the non-relativistic case [43l[T3] : 



1 



Dlk^ 



D. 



NR7,2 



Cp\{uj- csky + Df^^k^ (w + Csky + Df^^k^ 



where 



D. 



NR. 



,NR 



D 



,NR 



,(3.41) 

(3.42) 
(3.43) 



We see that relativistic effects appear only in the sound diffusion constant: 



D. = D. 



MR 



where 



MR _ 



(3.44) 
(3.45) 
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-0.1 -0.08-0.06-0.04-0.02 0.02 0.04 0.06 0.08 0.1 
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Figure 3.1: The spectral function for the Landau and minimal relativistic equa- 
tion. The solid and dashed lines respectively denote the Landau and minimal 
cases. The parameters are k = 0.1[l/fm], /i^ = 200[MeV], = 200[MeV], 
ij/i'iT'cSc) — C/{i^cSc) = 0.3 and \Tc/ (ncSc) = 0.6. Relativistic effects does not 
appear in the thermal mode but enhance the sound modes. 



First, the longitudinal kinetic viscosity is expressed in terms of the enthalpy 
density he in the relativistic case in place of the mass density see Eos. (13.341) 
and p.43p . We call this modification the minimal-relativistic (MR) effect. 

Next, the other is a genuine relativistic effect 5Ds which is absent in the non- 
relativistic case. This part comes from the dissipative term of Eq. (j3.10l) . which 
represent relativistic effects, and vanishes if we take the light speed c — ^ oo. 

To see the relativistic effects 5Ds quantitatively, we now determine thermo- 
dynamic quantities by the equation of state(EoS) of massless classical ideal gas, 
e = 3P = 3nT. Then, we have Cp = A, — 'i, ap — 1/Tc, Cg = •\/l/3 and the 
entropy density, Sc = — ficn-c/Tc. 

Figure 13.11 shows the spectral function, Eq. p.38|) . and the minimal rela- 
tivistic case with the above thermodynamic quantities. The parameter set 
is given as the following; k = G.l[l/fm], ^ic = 200[MeV], = 200[MeV], 
rj/sc — C/sc — 0.3 and XTc/sc ~ 0-6. 

As is expected, Fig l3. II shows that the sound mode is enhanced by the rela- 
tivistic effects, while the thermal mode is the same as in the minimal case. 

Now, we also compare our results to Eq. (|2.66p in which the coupling among 
Sn, Se and 6J is neglected. The spectrum (j2.66p has the only one slow mode, 
which is diffusive. In contrast, our result Ea. (|3.38p has the two modes: the 



and 




(3.46) 
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thermal and sound modes. In addition, we have the another slow mode in the 
transverse component Eq. p.30p : the viscous mode. 

Namely, the actual slow modes is three. In Sec l3.3[ we shall study the relative 
importance of them near the QCD critical point. 



3.2 For Israel-Stewart equation 

For the Israel-Stewart equation in the particle frame, the dissipative terms are 
T^"" = -nA^"' + + q^u^" + TT^"" , (3.47) 

and = {). Here 



(3.48) 



u^dpUv 



HP 



PiuPdpq, - aoa.n + aidpT:Pj , (3.49) 
2ryA^'''"^(apM, - hu'dr^pp - a^dpq^), (3.50) 



with u'^gp = 0, TT^"^ = it"!', u'^TT^t, = and tt^^ = 0. Here, ^o, Pi and /^a 
are the relaxation time of the bulk viscous, the heat flux and the shear viscous, 
respectively, ao (ai) is the coupling of the bulk viscose and the heat flux (the 
shear viscose and the heat flux). A^"'"^ is the projector defined by 



l^pvpa _ \ 

2 



-A'^^A'"^ 

3 



(3.51) 



Here, an important point is that the relaxation time corresponds to a correlation 
time of the memory function, Eq. (|2.9p . Namely, the Israel-Stewart equation has 
the memory effect, but such effect is irrelevant for the slow dynamics, as we 
shall see in the following. 

Applying the similar procedure as in the Landau equation, we have the 
spectral function of the density fiuctuation: 



S'n„(fe,Cj) 



(((5n(fc,t = 0))2) 



= 1- 



(w - CskY -t- £|2fc4 



2//3oC 



L^' + l/(/3oC)' 
2h,/[{f3ih, - l)XTo] 



{uj + cskf + Dlk^ 

+ 1/(2/327^)2 ^ ^2 + hl/[{(3^K - 1)AT, 



(3.52) 



We leave the detailed derivation to [26] . Here, we assumed the relaxation time 
satisfies the inequality, /3i > 1/hc- If the relaxation time does not satisfy the in- 
equality, we have a pathological behavior; the spectrum becomes negative. The 
Israel-Stewart equation in particle frame takes over the pathological behavior 
of the Eckart equation, in which the fluctuation does not relax; see the detailed 
discussion in [551. 



32 



3.3. TENDENCY AROUND THE QCD CRITICAL POINT 



Now, the spectral function apparently has six peaks including the conven- 
tional three peaks, but the new three Lorentzian functions should vanish in the 
long- wavelength limit fc — > 0, because the strength of these is of the second order 
in k. Therefore, the Israel-Stewart equation gives the same result for the hydro- 
dynamic spectrum as the Landau equation does in the long- wavelength limit. 
Namely, the relaxation times does not affect the result in the long-wavelength 
region. This result implies that the causality problem occurs only in the short- 
wavelength region. We note that such illegal component would be ruled out by 
the cutoff, as mentioned in Sec l2.1.3l 



3.3 Tendency around the QCD critical point 

Here, we include a part of information on the critical point in the hydrodynamic 
spectrum by static scaling laws. By this study, we shall find the tendency of the 
slow variables around the QCD critical point; namely, which fluctuations are 
enhanced near the critical point. We note that not all fluctuations are enhanced 
near the critical point. 

Now, we use the static scaling laws: 

-r/', (3.53) 
Cp (3.54) 

where ^ is the correlation length, which diverges at the critical point, a, v and 
7 are usual static critical exponents. The QCD critical point belongs to the 
static universality class of the 3d Ising model, Z{2). Then, we have the critical 
exponents 

a~0.1, (3.55) 
V - 0.6, (3.56) 
7^1.2, (3.57) 

and the critical behaviors of the specific heats are 

(3.58) 

Cp (3.59) 
Then, we have the critical behavior of the spectrum, Eq. (|3.38p . 



Snn{k,uj) I {{Sn{k,t = 0)y) = 



Cp/oj^ + D'^k^ 



D^k'^ D^k 



2 



Cp V (W - Cskf + L>2fc4 {uj + c,fc)2 + i:>2fc4 

2Dtk'^ 



(3.60) 
(3.61) 
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Figure 3.2: A schematic figure for the density spectrum near the critical point. 
The sound mode disappears, whereas the thermal mode is enhanced. 



because c„/cp behaves as and vanishes at the critical point. We see that 

the sound mode disappears, whereas the thermal mode survives. 
Moreover, the thermal diffusion constant behaves as 



A 



ncC, 



0. 



(3.62) 



Then, the surviving thermal mode is enhanced near the critical point. This is 
the critical-slowing down. 

Meanwhile, the transverse spectrum, Eq. (|3.30p . has no critical behavior. 
Then, the viscous mode is not enhanced nor suppressed near the critical point. 

Now, we see that relative importance of the three hydrodynamic modes. 

1. The thermal mode is the most relevant mode because it survives and is 
largely enhanced near the critical point. 

2. The viscous mode is the second relevant mode. The mode is not enhanced 
nor suppressed. 

3. The sound mode is suppressed and thus not so important. For the mini- 
mal critical dynamics, we can neglect the sound mode. However, for the 
renormalization of the bulk viscosity in the next chapter, we must consider 
the sound mode. An effect on the thermal mode by the sound mode is 
negligible, whereas a counter effect is large. 

From the above, we can know the tendency of the slow variables near the 
critical point. We first see that the density fluctuation, 5n, is enhanced by the 
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thermal mode, see Eqs. p.Gip and (I3.62[) . and Fig. 13.21 In addition, the energy 
density is also enhanced, because it also couples to the thermal mode. 

In contrast, we see that the momentum SJ is not enhance41, because the 
transverse and longitudinal components of that couple to only the viscous and 
sound modes, respectively. These modes are not enhanced near the critical 
point. 

Here, the critical behavior of the transport coefficients is not taken into 
account. However, we note that the similar analysis on the liquid-gas critical 
point, as in this section, qualitatively gives a good description [13]. 



^ We note that the momentum density is proportional to the fluid-velocity fluctuation: 
5 J = hcSv. 
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Chapter 4 



Nonlinear dynamics of the 
hydrodynamic modes by 
the dynamic RG 

Here, we study the nonlinear dynamics in the critical regime by the dynamics 
RG. We note that the nonlinear effects, which is not included in the usual 
hydrodynamics, arise in the critical regime. Thus, the analysis in the previous 
chapter is valid only in the hydrodynamic regime. 

In this chapter, we first construct the nonlinear Langevin equation to de- 
scribe the nonlinear dynamics in the critical regime. Next, we include the non- 
linear interaction in the transport coefficients by the dynamic RG. Thus, the 
transport coefficients are renormalized and diverge at the critical point. 

4.1 The nonlinear Langevin equation for the QCD 
critical point 

Now, let us construct the nonlinear Langevin equation. An important point on 
the construction is that Sn and Se are enhanced, while SJ is not near the QCD 
critical point, as mentioned in Sec l3.3l Then, we can neglect the nonlinearity of 
the momentum in the following construction. 

4.1.1 Thermodynamic potential for the slow variables 

First, we construct the thermodynamic potential H{5n, Se, SJ). 

Because the momentum density fluctuation is not enhanced near the QCD 
critical point, we now assume the potential for the momentum as the Gaussian 
form. Then, we have H(Sn, Se, J) — iJ„e((5n, Se) + Hj{SJ), with 

Hj{5J) ^ ^SJ\ (4.1) 
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In contrast to SJ, 6n and Se are enhanced near the QCD critical point, the ther- 
modynamic potential Hne{Sn,Se) should contain higher order terms of them. 

Now, the important point is that _ff„e((5n, (5e) is the quantity to determine 
the static property of the system and the QCD critical point belongs to the 
same static universality class as the 3d Ising class, namely, Z2. Therefore, we 
may construct Hne{Sn, Se) with the thermodynamic potential for the 3d Ising 
system [33], which reads 



l3Hisingiil^,m) = dr 



+joip'^m, H — 7=r™^ — hip — Tin 
2 Co 



(4.2) 



Here, ip and m are the spin density and the exchange energy density, respectively, 
ro, Kq^ uq, 7o and Cq denote the static parameters, while h and t the applied 
magnetic field and the reduced temperature, respectively. Then, we assume the 
thermodynamic potential as 

H{Sn, 5e, 5 J) = i/isinglV-, ™) + ^-^J^, (4.3) 

if we have the mapping between [tp, m) and (Jn, 5e). 

The general mapping relation between a grand canonical ensemble in Z2 
and the 3d Ising system is known in condensed matter physics [23], which are 
summarized as follows. First, we assume the following linear relation between 
the deviations of the intensive variables from the critical points, 

5h^ aid{fi/T) + a2ST/Tc, (4.4) 
6t^ /?i(5(^/T) + /32(5T/Te, (4.5) 

where ai , 02 , /3i and /32 are constants and assumed to be regular at the critical 
point. We note that ai, a2, /3i and /32 need not to be determined for the 
critical divergence, because those have no singularities at the critical point. 
Although one could use Eqs. (|4.4p and (|4.5[) for the mapping, it turns out to 
be inconvenient for a Langevin equation. To translate these relations to more 
convenient ones, we assume the following relation 

tPSh + mSr =::T-'^ST5e + S{fi/T)Sn, (4.6) 

which is actually derived by considering a change of the microscopic distributions 
by small deviations of the intensive variables in both systems. From the relations 
Eqs. (|4.4p - (|4.6l) . we arrive at the convenient mapping relation as follows, 

6n — aiip + j3im, (4-7) 
T-He^ a2^ + fi2m. (4.8) 



^ Recall that the static scaling laws are expressed by the deviations of the intensive variables 
from those at the critical point. 
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With this mapping, Eq. (j4.3|) now gives the thermodynamic potential for the 
QCD critical point. We note that we only map the static quantities, although 
the dynamic ones are studied. For later uses, we introduce fluctuations of the 
intensive variables as 

ST^T^-^, (4.9) 

This relation comes from the fact that, in the grand canonical distribution 
Pgra oc exp[(l/r)e+(/i/T)n], e and n are respectively conjugate to 1/T and /i/T 
[24] . We also introduce the fluid velocity fluctuation as in the non-relativistic 
case: 

r TT 

Sv^-. (4.11) 



We note that the static parameters in Eq. (|4.2p has the ultraviolet cutoff 
dependence in the region < A. Let us write the static parameters as ^(A), 
K{A), u{A), 7(A) and C(A) to make their A dependence, explicitly. These 
variables have the following asymptotic behaviors [27j Ell [33] : 



r(A) -A^-", (4.12) 

KiA) -A-", (4.13) 

u{A) ^A'-^"^, (4.14) 

7(A) ^ a(^+o'/^)/^-v^ (4.15) 

C(A) -A-"/^ (4.16) 

where e = 4 — c? with d being the space dimension, while a, v and t] are the 
usual static critical exponents. Noting that r\ is of order and very small, we 
neglect 77 and set Kq — 1, hereafter. 



4.1.2 Streaming terms and bare kinetic coefficients 

Here, we construct the streaming terms, u„, Ve and v j. We can nicely determine 
the first two terms from the continuity equations, because 8n and 8e are the 
conserved densities. From the continuity equations, we can write w„ and as 
divergences of reversible currents, which read 

in = (4.17) 
(e + P)72<5^;, (4.18) 

with j„ and jg being the reversible currents of the number and energy den- 
sity, respectively. Here, 7 is the Lorentz factor of the fluid-velocity fluctuation, 
n ^ rif. ^ bn and e = Cc -l- be. As the reference frame, we chose the rest frame 
of the equilibrium state, and then the back ground fluid velocity vanishes. Fur- 
thermore, We may set 7 ~ 1, because the fluid- velocity fluctuation is given by 
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Sv — ^5 J that is not enhanced. Therefore, we write the streaming terms, Vn 
and We, as 

Vn = -V • {n6v), (4.19) 
V, ^-V ■{{e + Pc)5v), (4.20) 

where we neglect the pressure fluctuation because it is not enhanced near the 
critical point |26j . 

Now, we note that the determination of vj is not simple. Although the 
continuity equation tells us that vj is the divergence of the reversible-stress 
tensor, the determination of the reversible-stress tensor is not trivial. However, 
we can determine it from the potential condition, which is a general condition 
for the streaming terms [24] • The potential (or divergence) condition [24l [13] 
rGcLcis 

j—n.e^J j—n,e,J •' 

We remark that this condition can be derived from Eg. (12.251) . In a continuum 
system, the right-hand side of Eq. (|4.21|) vanishes in general |24j . 
Thus, the potential condition is reduced to 

jdr Y: -^(^)^-O' (4.22) 

where vj is only the unknown quantity because we have already determined u„, 
Ve and H{Sn,Se,SJ). Using Eqs. ([ITT|) . ([il^ . (IT^ and g^l]), we obtain 



/ 



dr 



nV— + {e + Pc)V— -I- vj 
on oe 



(35v = 0. (4.23) 



Because this condition should be satisfied for an arbitrary fluid-velocity fluctu- 
ation, we have 

vj^-nV—-{e + P,)V—. (4.24) 
on oe 

Next, let us determine the kinetic coefficients from the relativistic hydrody- 
namic equation, Eqs. ([S21)-(|3IZ|)- From Eqs. ([32]), ([22]), ([iTTU)) and (|iTT|) . we 
can read the kinetic coefficients Lj^ for small 5v as 

nr-Tr^ ^ 



Lnn =-^o(^-pj (4.25) 
Vjj - -T,[rjQ5,,d,dj + (Co + (1 - 2/d)?7o)a.a,], (4.26) 

and that the other coefficients are zero. Here, d is the space dimension and the 
transport coefficients, Xo,rio and Co, are bare ones. 
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Now, we have determined all the terms, and then can write the nonlinear 
Langevin equation for the QCD critical point as 

-V.M.)-L„„« + ^„, (4.27) 
ot on 

^= -S7-iie + P,)Sv), (4.28) 

d6J ^SH , ^,^SH ^ S((3H) ^ 

— = -nV—-{e + P,)^—-Lj.j-^^ + e.j, 4.29 
ot on be oJ 

where 0„ and Oj are the noise terms and satisfy the fluctuation-dissipation 
relations 

{e^{r,t)e^{r\t))= -2Ao(^^) VH{r - r )5{t - t (4.30) 

{e:j{r,t)e'j(r\t)) = -2r,[r,o^'^V2 + {Co + (1 - 2/d)r^^}d'&] 

xd{r -r')6{t-t). (4.31) 

Here, an important point is that the first and second terms in Eq. (|4.29p denote 
the nonlinear effects that are absent in the usual hydrodynamics. These terms 
represent the time-reversible forces acting on the fluid, and the force in the usual 
hydrodynamics is only the pressure gradient, —VP. Thus, the nonlinear part 
in those terms reflect the softening of the thermodynamic potential and is not 
included in the hydrodynamics. 

Let us now write the transport coefficients as A(A), r]{A) and C(A) to make 
their cutoff dependence in the critical region. The critical behaviors of the 
transport coefficients are determined from their asymptotic behaviors near the 
relevant fixed point as A is lowered. 

Now, we compare the Langevin equation, Eqs. (|4.27p - (|4.29p . with that for 
the liquid-gas critical point [37] 

d6n 
~dt 

^ = -V • ((e + P,)5v) + AoT.V^^ + 0^, (4.33) 
at oe 

^ . _„vf!s - (. + Pjvfi - L.. . M + (4,34) 

ot On oe oJp 

where SJp = PcSv, p and H^^ are the non-relativistic momentum density, the 
mass density and the thermodynamic potential for liquid-gas critical point, re- 
spectively: 

1 2 

H\g{Sn,Se,6Jp) = _ffising('0, m) -I- - — 6 J p. (4.35) 

We see that the streaming terms have the same forms but the dissipative ones are 
different between the relativistic and non-relativistic cases. The difference also 
appears the relation between the momentum and the fluid- velocity fluctuation. 
Therefore, one may naturally expect some novel characteristics in the relativistic 
case that is absent in the non-relativistic case i27i . 



-V-{n5v), (4.32) 
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4.2 The transport coefficients by dynamic RG 

Here, we study the transport coefficients by the dynamic RG. A detailed deriva- 
tion of the RG equations is given in the Appendixes. 

First, we rewrite Eqs. (|4.27l) - (|4.29p as the equation in terms of {ip,m), 
to conform (6n,Se) to (ipjrn). Noting that we can set 0:2 = in the mapping 
relations, Eqs. (|4.7|) and (|4.8|) . without loss of generality (27!, have 

^= -C^V.5J-a^^Kj-^^~h-'V-{^jSJ)+a^'9r„ (4.36) 



dm 
'dt 



I3:^'V ■ SJ - h-'V ■ (mSJ), (4.37) 



d{SJ) _ SH ^-,.^6H 

SH SH 

_^V— - mV— - {T,K)~^Ljj ■ SJ + Oj, (4.38) 

dip dm 

with = a^^{nch~^ — /3i/3^^) and Cj = a^^{nc — Pihc)- Here, we note that 
we could rewrite the potential, Ea. (|4.2l) . as that in terms of {6n, Se) to conform 
the variables; the choice is a matter of preference. 

In the dynamic RG transformation, we average over the short wavelength 
components in the shell, A — (5A < fc < A, for the Langevin equation. For this 
task, we must perturbatively solve the equation about them, by rewriting it 
as a self-consistent equation [13 . Although an explicit derivation of the self- 
consistent equation for the QCD critical point is first made in this thesis, we 
leave the details of the derivation to Appendix O because the general procedure 
of the derivation is standard and given in the textbook [13]- Here, we provide 
only a few basic equations of the dynamic RG for the QCD critical point. 

Now, as is shown in Appendix [Cl Eqs. (|4.36p - (|4.38l) can be reduced to the 
following form; 

/^(fe,c.)\ /^°(fc,w) 
m(fe,w) = ( I +G°{k,io)V{k,Lu), (4.39) 



and 



SJ^{k,uj) = Sj"^{k,uj) + G"^V^^^{k,Lo), (4.40) 



where 5J||(fe) = k ■ 6J{k) and 6J±{k) = 5J{k) — 5J\\{k,) are the longitudinal 
and transverse components of the momentum. Here, G*^ and G\ are the bare 
propagators, which are given by Eqs. (|C.50|) and (IC.51|) - (IC.53I) . whereas 
V and Y the nonlinear couplings, coming from the streaming terms and 

given by Eqs. (Ic:29]) - (IClM)) and Also, , SJ^^ and 6J°j_ are the bare 

variables, which are the solutions without the nonlinear terms. Iterating the self- 
consistent equations (14.391) and (j4.40p . we can obtain a perturbative expansion 
of the nonlinear couplings and have a coarse-grained Langevin equation. 
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Now, we note that the variables, i^^Ji and J||, are respectively correspond to 
the thermal, viscous and sound modesG (see the propagators (IC.51|) - (IC.53|) .). 
Thus, the first and third rows of Eq. (|4.39l) respectively denote the equations 
of motion for the thermal and sound modes, while Eq. (j4.40p for the viscous 
mode. We stress that the sound mode is neglected in the model H, although it 
is essential for the renormalization of the bulk viscosity. 

Here, we make the coarse graining to the second order in the nonlinear 
couplings, V and V±^^ (see Fig. ID. 31 for an example.). Inspecting the coarse- 
grained equation for (see Eq. (ID.5|) for the detail), we have the RG equation 
for the thermal conductivity: 

-A^^ J/(A)A(A), (4.41) 

/(A) = TcK4/{D^ri{A)X{A)A'^), K4 is the surface area of a unit sphere in 4 
dimensions divided by (27r)'*, = {ncT,./ aihcY' . Here, we have introduced 
/(A) for convenience sake. Similarly, from the coarse-grained equations for 
and J||, we obtain the RG equations for the shear and bulk viscosities 

_^^C(A) ^ ^^2/^^^-l.^w-.-4^ (4 

oh 

where 7(A) is a static parameter in the thermodynamic potential (see Eqs. (|4.2p 
and (|4.15p ). and A = hlK^/ {I32D^). Furthermore, differentiating /(A) about 
A, we also have the RG equation for it: 

-A2M.;,^)(,_g;,^,)^ ,«4, 

Now, we note that Eqs. (|4.41l) . (|4.42l) and (|4.44l) are identical to those for 
the hquid-gas critical point except for unimportant constants in f (A)[24l [27]. 
Equation (|4.43p is also equivalent to the RG equation of the bulk viscosity for 
the liquid-gas critical point in the limit w — [Ml [57] . Therefore, arguments 
about the RG equations and results from those are the same as for the liquid- 
gas critical point. Then, we provide only essential arguments and results in the 
following part, and leave the detail to [^ [5^ [511 [?7] . 

Now, we identify the relevant-fixed point as the following [24l[23. Because, 
at a fixed point, parameters are invariant about the RG transformation, we set 
the left-hand side of Eq. (|4.44|) as 0. Then, as a fixed-point value of /(A) which 
is denoted by /*, we have /* = and /* — (24/19)e. Therefore, we have the 
two fixed points, and the relevant one is specified by /* = (24/19)e. Although 
the relevant point seems to be absent in Eqs. (|4.4ip . (|4.42p and (|4.43p . the 
reason is due to the simplified RG transformation as mentioned in the earlier 
section, and this is just a apparent problem [311 [SI] ■ 



^Although m would be a linear combination of the thermal and sound modes, we need not 
to consider m for a following analysis. 
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Substituting /* = (24/19)e into Eqs. ^ZH, and we have the 

asymptotic behaviors near the relevant-fixed point: 

(4.45) 
(4.46) 
(4.47) 

Here, in the derivation of Eq. (|4.47p . we have used the asymptotic behavior of 
7(A), Eq. ()4.15p . Decreasing the cutoff to the region A <^ we can replace 
A with in the asymptotic behaviors [24j l3T|: 



A(A)^ 


^ A- 




ry(A)^ 


^ A- 




C(A)^ 


^ A~ 


-(4- 



Ar - (4.48) 
77R ^ (4.49) 



Cr^^'"^^~^- (4.50) 



In three dimensions, we find 



Ar ^ (4.51) 

^R - (4.52) 

Cr - C'-'. (4.53) 

We can also read the dynamic critical exponents from Eqs. (|4.48p - (|4.50p . A 
dynamic critical exponent, denoted by z, generally parametrizes the decay rate 

r(fc) at the wavenumber k — as r(^~^) ~ As shown in Appendix [Cl 
the decay rates for the three modes at k are given by 

^ thermal 

(fc) = Xuk^{rR + k^)D^, (4.54) 

rviscous(fc) = ■quk'^h^^, (4.55) 

r.ou„d(fc)= (CR + 2(l-l/d)r7R)A:2/i-i. (4.56) 

Thus, we find the dynamic critical exponents as 

^thermal = 4 - ^6, (4.57) 

'^viscous — 2 777^' (4.58) 



19 

19' 



•-sound = - ( 2 - - - ) . (4.59) 



In three dimensions, those are 



^thermal 3, (4.60) 
■^viscous 2, (4.61) 

2sound 0.8. (4.62) 
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We see that the thermal and viscous modes exhibit critical slowing down, while 
the sound mode critical speeding up. 

Why do not the relativistic effects appear in the RG equations? The reason 
is that the nonlinear terms in the dissipative terms generally renormalize only 
static parameters, up to order [31|, I13j. Furthermore, the difference in the 
relation between the momentum and the fluid velocity is only unimportant 
constants, i.e., the enthalpy density h and the mass density p. Then, the RG 
equations are essentially the same as for the non-relativistic case. 
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Chapter 5 

Summary and Concluding 
remarks 

We studied the critical dynamics near the QCD critical point by the linearized 
relativistic hydrodynamics and the dynamic renormalization group (RG). 

First, we studied the linear dynamics of the baryon number and the energy- 
momentum in the hydrodynamic regime, fc -C C^^, by the relativistic hydrody- 
namics. We showed that the actual slow modes are three: the thermal, viscous 
and sound modes. Furthermore, we found that, near the critical point, the 
thermal mode is enhanced and the most relevant. In contrast, the sound is 
suppressed and negligible for a minimal critical dynamics. 

By this study. We also found that the Landau equation, which is believed to 
be an acausal hydrodynamic equation, has no problem to describe the slow dy- 
namics. We also showed that the Israel- Stewart equation, which has relaxation 
time and is causal, gives the same slow dynamics as the Landau equation gives. 
This result implies that the causality problem occurs only the short-wavelength 
region. We note that the short-wavelength region would be out of applicable 
scope for the hydrodynamics and actually ruled out by the cutoff. We also stress 
that the relaxation time is important only for a rapid motion. 

Next, we studied the nonlinear dynamics in the critical regime, <C k, 
by the dynamic RG. For this purpose, we constructed the nonlinear Langevin 
equation near the critical point for the first time. Our construction is based 
on the generalized Langevin theory, by Mori [351 123] , and the relativistic hy- 
drodynamics; instead of a naive construction method |13| . we determined the 
streaming terms by the relativistic hydrodynamics and the potential condition, 
which gives a constraint to these terms. The resulting equation is given by Eqs. 
(|4.27|) - (|4.29|) . Although there are some attempts to make a one-to-one mapping 
between the QCD critical point and the Ising critical point [iQlISS], we showed 
that it is not necessary to specify such the mapping for the critical exponents, 
as for the liquid-gas critical point [24] . 

We showed that the bulk viscosity and the thermal conductivity strongly 
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diverge at the QCD critical point. Also, we found that the thermal and vis- 
cous modes exhibit critical slowing down with the dynamic critical exponents 
•^thermal ^ 3 and Zviscous ^ 2, respectively. In contrast, the sound mode critical- 
speeding up with the negative exponent Zsound ~ —0.8. We stress that the 
earlier studies [HI [21] treat the bare transport coefficients and does not include 
the macroscopic nonlinear interaction. 

We note that the bulk viscosity and the thermal conductivity are usually 
neglected in heavy ion physics, however they become much more important than 
the shear viscosity near the QCD critical point. Furthermore, the description 
for the created matter as a perfect fluid is not valid near the QCD critical point 
by the strong divergence of the bulk viscosity. 

As the argument about the dynamic universality class [20|[T4], we showed, 
from an explicit calculation, that the QCD critical point has the same criti- 
cal behaviors as the liquid-gas critical point has. The argument assumes the 
insignificance of the relativity for the critical dynamics by the slowness of the 
diffusion processes. However, we showed that the genuine reason for the in- 
significance originates from the small fluctuation of the momentum density; the 
critical dynamics is essentially governed by the streaming terms, which are mod- 
ifled by the relativistic effect through only a Lorentz factor of the fluid velocity 
fluctuation. However, the fluid- velocity fluctuation, which is proportional to the 
momentum, is not enhanced near the critical point. Thus, the relativistic effects 
do not affect the critical dynamics near the critical point. We stress that the 
sound mode exhibit critical speeding up, and then the sound diffusion is fast 
near the critical point. Therefore, the basis of the conjecture would be true for 
the thermal and viscous modes, but not for the sound mode. We also note that 
the model H [22], which is the minimal-dynamic model for the dynamics near 
the liquid-gas critical point, can not describe the critical behavior of the bulk 
viscosity because it does not contain the sound mode. 

We note that our Langevin equation must satisfy usual fluctuation-dissipation 
relations, Eqs. (|4.30p and (|4.31|) . for the consistency with the linearized Lan- 
dau equation 0, although a relativistic Brownian motion seems not to satisfy 
the usual relations [33]. Moreover, our Langevin equation seems to violate the 
causality, because the dissipative terms are determined from the Landau equa- 
tion. However, the Israel- Stewart equation, in which the causality problem is 
formally resolved, gives the same result as the Landau equation gives in long- 
wavelength region, as shown in Sec l3.2l Therefore, our determination from the 
Landau equation must suffice. 

Also, we note a frame dependence of the results in Chap. H] As a hydrody- 
namic equation, we used only the equation in the energy frame. Does the results 
change if an equation in the particle frame is used? Although the frame depen- 
dence can appear in only dissipative terms, the critical dynamics is essentially 
determined by the streaming terms. Therefore, the results would not change for 
the particle frame, if an equation in the frame is correct. However, in practice. 



^If our nonlinear Langevin equation is linearized, the linearized equation must give the 
same result as the Landau equation gives. 
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the Eckart equation has a pathological behavior [TB]. Namely, fluctuations do 
not relax, and therefore we cannot use the Eckart equation. 

Furthermore, we note Lorentz covariance. Our Langevin equation is not 
Lorentz covariant, but it would not be a problem. The reason is the following. 
Here, we consider the fluctuations in the background medium. In such situation, 
the Lorentz transformation boosts the fluctuations but not the medium. Then, 
after the Lorentz boost, we have the boosted fluctuation and the medium that 
still rests. Namely, the boosted system differs from that before the boost. The 
covariance means that, if we see the same system from different reference frames, 
we have the same physics. Thus, our Langevin equation would be no problem. 
Actually, an equation of the Brownian motion is not Galilei covariant by the 
same reason, but it is no problem. 

Recently, some authors suggest the existence of other critical points in higher 
density region of the QCD phase diagram where the color superconductivity 
is taken into account |42j [41]. It would be interesting to study the critical 
dynamics near such a new QCD critical point using the dynamic RG theory, as 
an extension of the present work. For this purpose, however, we must firstly 
specify the soft modes and construct the nonlinear Langevin equation. If the 
soft modes are different from the conserved densities, which is the case when 
the diquark fluctuations are relevant [42l |37l |44] , the construction based on the 
relativistic hydrodynamics done in the present work does not work, and we must 
directly recourse to Eq. ()2.47p to identify the streaming terms. 
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Appendix A 

Generalized nonlinear 
Langevin equation 



Here, we give the detailed derivation of the nonlinear Langevin equation (|2.24l) . 
As mentioned in 12.1.21 we can obtain the nonlinear Langevin equation about 
Aj from the linear Langevin equation about g{A,a) [MllS]- 

Thus, let us first derive the linear equation about g{A,a). We note that the 
operator identity Eq. (|2.5|) is valid even for Vg. Multiplying Ea. (j2.5|) by g{A, a), 
we have the linear Langevin equation for g(A(t), a) = exp[iCt\g(A, a): 



^^g{A{t),a)^ J dain^^,g{A{t),a)-j dt J da^^^>it)g{A{t-t),a)+Fa{t). 

(A.l) 

Here, 

in^^, ={{zLS{A-amA-a))/P,,{a), (A.2) 



where Vj{a) is 



The noise Fa{t) is 



vjia) = {AjSiA - a))/P,^{a) = (i,;a). (A.4) 



Fa{t) = cMQg^Lt]QgiL5{A- a) (A.5) 

= -Y, ^[c^V [QaiLt]{QgAu)5{A - a)] (A.6) 
k 

k 
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where we introduced 



U{t) ==exp[QgiLi], 
ej{t) = U{t)QgiLAj{Q). 

With the noise, the memory function {t) is written as 
^aa'{i)^{Fa{t)FAQ))/P,,{a). 



(A.8) 
(A.9) 



(A.IO) 



If we multiply Eq. ljA.ll) by Oj and integrate over a, we can have an exact relation 
corresponding to Ea. (|2.6p . 

However, let us now make a Markov approximation because the exact ex- 
pression is inconvenient for our purpose. Again, if the time scale of Aj is much 
larger than that of we can make the approximations: 



d 



and 



dt J da^^^>{t )g{A{t-t),a)r 
Introducing a bare transport coefficients as, 



da 



(A.11) 



g{A{t),a ). 

(A.12) 



Lik{a) 



dt{9,{t)ek{0);a). 



we can rewrite the memory function as 

1 



dt^aa'W 



Pcq{a') Jq 



dt{Fa{t)F^-{0)), 



(A.13) 



(A.14) 



7^ E J-^iL^k{a)P,,{a)S{a - a )]. (A.15) 



Peg (a') ^ da, da 



After these approximations, we can obtain the nonlinear Langevin equation 
^^n^ from Eq.dXT]). 
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Appendix B 



Detailed derivation in the 
chapter S 



Here, we gives the detailed derivation of Eq. p.38p . 

Let us first calculate the matrix elements, {A^^)ii and (A^^)i3, in Ea. (l3.35p . 
They are given by the simple formula (A^^)ii = gj^(>i22^33 — ^23^32) and 
(A-i)i3 = g-L^(Ai2^23 - A13A22). Here, det A reads 



det A 



TlcCy 



z'k' ( ^ A + 1^1+ - 2DtclapT, 

Cy fir 



(B.l) 



where denotes the higher order terms in k. We are interested in the low- 
wavenumber region. Then, det A can be nicely factorized to second order in 

""""■■-(z + Dtk'^){z + Dsk^ + icsk){z + D^k^ - ic.fc), (B.2) 



where 



hrTr 



Ds = - 
2 



(B.3) 



Then, we can write the Fourier-Laplace coefficient of the density fiuctuation 
to second order in fc, 

5n{k, z) (2 + Dsk"^ + icsk){z + D^P - ic^k) + zk'^^^f^ - fc^^ 
Sn{k,0) {z + Afc2)(z + + icsk){z + D,k'^ - icsk) ' 

Performing the inverse Laplace transformation 

dn{k,t)/6n{k,0) ^ — dze^'^Snik, z)/6n{k,0), (B.4) 

27''i Js-ioo 
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we obtain the dynamical density fluctuation at t 



6n{k,t)/5n{k,0) 



— cos{cskt)e' 



(B.5) 



Here, we have retained only the terms in the amplitudes to zeroth order in 
k. Because Ea. (|B.5l) is the density fluctuation in a stationary process, we can 
replace the time t hy \t\. Therefore the Fourier transformation of Eq. ljB.sp is 



Sn{k,uj)/Sn{k,0) = 1 - 



Cp [(a;-c,fc)2 + L>2/c4 {u + Cskf + D^,k^ 
Thus, we finally obtain the spectral function of the density fluctuation 
Snn{k,uj) — {Sn{k,uj)dn{k,t — 0)) 



{{Sn{k,t = 0)f) 



1 



Cpjuj^ + D^k* 



Cp\{uj - Cskf + D^k^ {lj + Cskf + Dlk^ 



(B.6) 



(B.7) 
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Appendix C 

Rewriting the nonlinear 
Langevin equation as a 
self-consistent equation 



Here, we rewrite the Langevin equation, Eqs. ()4.36|) - ()4.38|) as a self-consistent 
equation. First, we make a Fourier transformation as the foUowing 



Then, we have 



2r~ 



-C^ik ■ SJ - Lnn 
-h-^ik- ( {^:{q)5J{k~q))+a~^t 

JqVt 



—iujm(k, uj) 
-iujSJ{k, u) 



^{q)ipik -q) + ^{q)m{k - q) 
dip dm 



-{T,h,)-'Ljj-SJ + e, 



(C.l) 



(C.2) 



-p^hk-dJ -K^ik- / {m{q)SJ{k-q)), (C.3) 
Jqn 

-Cjik^ ~ (i^^hcik^ 
dip dm 



(C.4) 



Note that the quantities with tilde in Eq. (jC.2l) - (jC.4p are Fourier transformed, 
like Eq. (jC.ip . and we have abbreviated the nonlinear terms such as 



lJj{q)SJ{k - q) 



qfl 



dn d'^q ~ 



27r {2tt' 



^^iq,n)SJ{k-q,u-n). 



(C.5) 
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AS A SELF-CONSISTENT EQUATION 



We now decompose Eq. (jC.4[) into the longitudinal and the transverse compo- 
nents: 



j5J\\ = 



-iuj5J\ 



dip dm 



-I / (k-q) 



^{q)ip{k -q) + ^{q)m{k - q) 
dip dm 



-{T,H,)-^k-Ljjik)-dJ + e\i, 
'SH 



-I I r^{k)-q 



dH 



{q)ijj{k - (?) + — {q)m{k ~ q) 
dip dm 



iTcK)-^V^{k) ■ Ljj{k) -dJ + Oj 



(C.6) 



(C.7) 



where we introduced a projection operator as 



{V±{k))ij = dij - k.kj/k^, 



(C.8) 



and 



5J\\{k) =k-6J{k), 
SJ_L{k) =V±{k)-6J{k), 
e\\ik) =k-e{k), 

e±{k) =vAk)-e{k). 



(C.9) 
(C.IO) 
(C.ll) 
(C.12) 



Because the streaming terms in Eqs. (|C.6I) and (|C.7p are too complicated for 
our purpose, let us retain only the terms that yield dominant contributions for 
the transport coefficients. We note that only such terms suffice for obtaining 
the critical exponents. From the relations 



d''r{^{r)iP{0)) 
dV(m(r)m(0)) ~ 



(C.13) 
(C.14) 



we expect ijj yields stronger singularity than m. Therefore, we only retain the 
term that are of the second order in ip. Namely, we reduce the streaming terms 
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as 



dip dm 
SH 



+i (fe ■ q) 

Jqn 



SH , 



^^iq)^ik ^ q) + —{q)m{k - q) 



+1^2 ^hck^o / tjjiq)ip{k - q) 

Jqfl 



Jqn 



Yt(i)^(^ - g) + ^{q)rh{k - q) 
dip dm 



iT,V^{k)- / qx^'{qmq)4>{k-q). 
Jqn 



(C.15) 



(C.16) 



where Xo ^(k) = + k^. Notice that we set Kq = 1, as mentioned in the text. 

Next, we consider the dissipative terms. The important point is that the 
nonhnear terms in dissipative terms generally renormalize only static parameters 
in a thermodynamic potential, up to second order in e [311 113] . Therefore, we 
can take into account nonlinear terms in the dissipative terms with the results of 
static RG, Eq. (|4.12l) - (|4.16l) . and effectively neglect it in the Langevin equation. 
Then, we reduce the Lnn5{(iH) / dip as 



(C.17) 



In contrast, the dissipative terms of 5 J are originally linear and then directly 
read 



k-Ljj{k)-SJ{k,uj) ^T,[Co + 2{l-l/d)Tjo]k^6Jn{k,u), (C.18) 
P^{k) ■ Ljj{k) ■ 6J{k,Lu) ^Tcriok^SJj_{k,uj). (C.19) 



Collecting the above results, we arrive at the reduced nonlinear Langevin 
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APPENDIX C. REWRITING THE NONLINEAR LANGEVIN EQUATION 
AS A SELF-CONSISTENT EQUATION 



equation: 



— iLU-il) = —ikC^5J\\ — ^ik ■ / 'ijj{q)SJ{k — q) 

Jqn 

—iujfh = —P^^ikdJ^ — h~^ik ■ / rh{q)6J{k — q), 

Jan 



(C.20) 
(G.21) 



—iuSJn = Tc 



~ikP2 ^^c7o / ^{q)ip{k - q) 



-iuj6J±= -inV±{k)- [ qXo^QMlMk - q) 

JqQ 



where 



-k'^mK ^sj± + e±, 



ncTc 



v\ =[Co + 2(l-l/%o]. 



(C.22) 
(G.23) 

(C.24) 
(C.25) 



This is the basic equation for the dynamics near the QGD critical point, which 
is first written down, and a main result of this paper. 

We can compactly rewrite the basic equation in a matrix form: 



where 



M(k,i^) 



M.(k,LS) I m{k,ijj) 
\5J||(fe,a;) 







V{k,uj) + e{k,uj), 



(C.26) 





-iu) 



ik/3^' 



ikxo\k)Cjn ikC^^p2hcTc -iuj + k^ul,h-\ 
e{k,u)) 



^a^^e{k,Lo)^ 


0u{k,uj) 



(V^^± {k,uj) + V^^w {k, uj) 
Vmm±{k,u) + V^^\l{k,u) I , 
V\\^^{k,u;) 



(C.27) 



(C.28) 



(G.29) 
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and 



= ~h^^ik- / ip{q)SJ_i_{k~q), 

= ~h-^i [ k-{k-q)/\k-q\ 
Jnn 



xi'iq)SJ\\{k-q), 



-h^ ^ik 



■m{q)SJ ±{k — q), 



ik / ■m{q)SJ\\{k~q), 
'ikTcf^2^hcjo '4'{q)i^{k ~ q) 

JqO, 



(C.30) 

(C.31) 
(C.32) 

(C.33) 

(C.34) 



Because Eq. (|C.23P is decoupled from the other equations at hnear level, we do 
not rewrite it as the matrix form. 

Next, we calculate the bare propagator G'^{k,uj) = j\4^^{k,uj). The in- 
verse matrix is given as the transposed cofactor matrix divided by det The 
determinant reads 

detX = i-iojf + {-iujfk^{XD^Xo\f^) + ^oK^) 
^iujk\C„^hcT, + Xo\k)C^CjT,) 
+k^XoXo^D^,Ca^h,T, 

-iojk^XoXQ^D^,{k)voK^. (C.35) 

Here, in the coefficient of —iujk'^, taking into account the behaviors after renor- 
malization 1271 [33l , which are 



0.2 



(C.36) 
(C.37) 



we neglect Xo ^{k)C^C,jTc by comparing with Cq ^hcTc . Then, we can factorize 
the determinant as 

detX - + Ao(fc)xo 



X {-iio + ikcs + -^i^oh^ k^) 
x{—iuj — ikcs + ^VQh^^k^), 
in the long-wavelength region. Here, we defined 

Ao(fe) = Aofc^-D^, 



(C.38) 



(C.39) 
(C.40) 
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APPENDIX C. REWRITING THE NONLINEAR LANGEVIN EQUATION 
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The diagonal components of the cofactor matrix m reads 
mil {—ico + ikcs + ^Poh~^ k^) 

x{—iuj — ikcs + ^^'o^c^k'^), 

m22 = (-i<^)^ - i^^k'^>^oXo{k)D^i/^Qh~'^ 
+k\o\k)C^CjT,, 
+k^\oXo{k)D^y„h-^, 
= + Ao(fc)xo^(fe)), 

and the off-diagonal components are given by 



(C.41) 



(C.42) 
(G.43) 



mi2 


= k\^\k)/3^'CjT,, 


(C.44) 


mi3 


= -kuxo\f')CjTc, 


(C.45) 


m2i 


= — fc^Cg ^hcC^f32Tc, 


(C.46) 


m23 


= -ikCo^hep2Tc{-ioj + e\oXo{k)D^), 


(C.47) 


msi 


= k'^CQ^hcC^p2Tc, 


(C.48) 




= -ikl3^\-iu> + Xo{k)xo\k)). 


(C.49) 



Here, we neglect the off-diagonal components because they would not yield 
dominant contributions to the transport coefficients. Then, we obtain the bare 
propagator as 



with 



G°{k,oj) = 

G%{k,uj) 
G\{k,uj) 



G%{k,u:) 





Gl,{k,Lo) 
G\{k,uj). 



-iu: + \{k)xo\ky 



1 



-iijj + ikcg + ^i^Qhc^k"^ 
1 



—iij — ikcs + ^f^hc ^k^ 



(C.50) 



(C.51) 



(C.52) 



Gmm{k,'jj) is not needed in later calculations. The bare propagator of 6J± is 
trivially given by 

Gl{k,uj)= . J ,,,_i . (C.53) 
We finally arrive at the equations of motion as the self-consistent form: 



■0(fc,WJ 

m(fc, w) 
\^SJn (fe, w) ; 



■ V3"(fc. 


^<5jO(fc,C.); 



+ G°(fe,u;)y(fe,w), 



(C.54) 
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and 

SJ^{k,Lo)^Sj]_{k,u)+GlV^^4k,u), (C.55) 

where 

^\k,uj) ^Gl{k,Lu)a^^0n{k,Lu), (C.56) 
dJ^^{k,uj) = G|j(fc,w)^||(fc,w), (C.57) 

Sj\{k,u) =G'i{k,uj)e^{k,u), (C.58) 

V^^4k,u;)^ -iT.VAk)- f qXo'iQ)HQ)i'ik-q). (C.59) 

Here, ■ip^{k,uj), SJ^^{k,uj) and Sj'^^ are the bare variables that are the solutions 
without the nonlinear terms. Iterating Eqs. (|C.54p and (jC.55|) . we can obtain 
perturbative expansions about nonlinear interactions V and V ±_^^. We note 
that the first and third rows of Eq. (jC.54|) are the equations of motion for the 
thermal and sound modes, respectively, while Eq. (jC.55|) is for the viscous mode. 
We also stress that Eqs. (|C3T|) - ([C33l) are the propagators of the thermal, sound 
and viscous modes, respectively. 

Now, we calculate the two body correlation of 'tlP{k, lo) and 5J^ (fc, w), which 
are needed in later calculations. 

(Vi°(fei,^i)Vi°(fc2,c^2)) = GS),(fei,L.i)G^(fc2,a.2)«r' 

X(^(fci,wi)^(fc2,^2)). (C.60) 

Using the fluctuation dissipation relation Eq. (|4.30p . we find 

{e{ki,L0i)e{k2, L02)) - 2alXo{k){27rf+^Sih + fc2), (C.61) 



and 



{^%k„.,)^'{k,,.,))= 

+ A^(fci)xo (fci) 
x{2^f+^5{ki+k2), (C.62) 



(Wl(fei,^i)5Jl(fe2,^2)) - 2^% .. ,-2 (^^(fci))^ 



where ^(fci + k2) = (5(fci + k2)5{uJi + 0-12)- By a similar calculation, we have 

2TrM{ki) 
ujI + ril{ki)hc 
x{2nf+^Siki+k2), (C.63) 

where r]o{k) — rj^k^ . For a later convenience, we define 

2Ao(fc) 

uj'' + xl{k)x^'\ky 
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: 





^ J I : SjV. G-o^ : - ^ - : - O - 



Figure D.l: Diagrams for the full and bare variables, the bare propagators and 
the bare correlations. 



Appendix D 

Renormalization of the 
thermal and viscous modes 



Here, we first derive the RG equations for the thermal conductivity and the shear 
viscosity. Now, we note that the sound mode is not a genuine-relevant mode but 
a secondary mode that is strongly affected by order-parameter fluctuations but 
yields only a negligible feedback for the order parameters [261 13S] • Then, we can 
neglect the sound mode for the minimal critical dynamics; however, the bulk 
viscosity is not renormalized in that case. Here, to first analyze the minimal 
dynamics, we neglect the secondary mode, which is renormalized in the next 
section. In that case, the equations of motion are 



and Eq. (IC.SSp . For a diagrammatic treatment , we denote the full and bare 
variables, the bare propagators and the bare correlation functions as Fig. ID. 11 
Then, we can represent the equations of motion (jD.ip and (IC.55|) as Fig. ID. 21 

For coarse gaining, we decompose the variables into the long- and short- 
wavelength components as 



(D.l) 



(D.2) 
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< - --< 

Figure D.2: Diagrams of the equations of motion for the thermal and viscous 
modes. The left and right hand side respectively denote Eqs. (jD.ip and (jC.SSp . 




Figure D.3: Diagrams for the coarse-grained equation of motion for ip. The 
letters, L and S, respectively denote the long- and short-wavelength components 
(see the text below Ea. (|D.4p ). 



with 

i;^{k,uj)=Q{A-SA-k)ij^{k,uj), (D.3) 
4>^{k,uj)= Q{k-A-SA)i^{k,uj), (D.4) 

where Q{x) is a step function; i.e., the wavenumber is decomposed into < 
k < A — 6A and A — 6A < k < A. Hereafter, quantities with the suffixes L 
and S are supposed to be decomposed as above. To average over the ip'^^ and 

~ OS 

SJ^ , we must solve the equation of motion about them. Here, we solve the 
equations to second order in the nonlinear interactions and average over ■0"^ 

~ OS 

and 6J j_ . Then, we find the coarse-grained equation of motion for ip, which is 
diagrammatically given by Fig. ID. 31 The last two terms in Fig. ID. 31 represent 
nonlinear interactions being of third order, and can be neglected. Furthermore, 
the fifth term vanishes due to the relation between the step and delta functions 
in the loop integral. Introducing the self energy E^^, which is graphically 
represented in Fig. ID. 41 we can write the coarse-grained equation for as 

+i}^{k,uj)Gf{k,io)^^^{k,uj). (D.5) 



64 



Figure D.4: Diagrams for the self energies. 



The self energy is 



{k-r{k-q)-k)xoiq) 



(D.6) 



/, -iuj + Ao(q)xo ^(9) + Voik- q)hc ^ ' 
where r]o{k) — ijok^. Solving Eq. (jD.5|) about ip^, we have 

+ [{G°^)-^-j:^^r'V^^^. (D.7) 

where we used Eq. (jC.56p . Introducing renormalized variables as 

(Gv,R)-i(fc,c^)= {Gf)-\k,u)-i:^^{k,uj), (D.8) 
iP'^^{k,u)^ GMk,Lj)a^'e„{k,u), (D.9) 

we can rewrite Ea. (jD.6p as the renormalized equation of motion: 

= i^°^{k,u) + GMI^,uj)V^^^. (D.IO) 

We now require that the renormalized propagator has the same form as the bare 
one: 

{G^n)-\k,uj) = ~iLj + \RD^ex-\k), (D.ll) 

where Ar is the renormalized thermal conductivity. That is, we require that the 
only transport coefficients are explicitly renormalized. The small correction for 
the thermal conductivity 5A = Ar — Aq reads 

5\ = - lim [{D^k^xo{k))~'^T.^^p{k,u:)\, 

k.UJ—¥0 



f {k-V{q)-k)xo{q) 



hcD^ Jg Xo{q)xo ^(.q) + Vo{q)hc ^ 
We approximate the denominator and the numerator as 



(D.12) 



Mq)Xo'iq) +Vo{q)h-^ - Vo{k)h-\ (D.13) 
Xo\q)^ ro + q^ ^q^, (D.14) 
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near the critical point |31) . Then, we find 



where dQd is the sohd angle in the space dimension d. Therefore, we obtain the 
RG equation for the thermal conductivity: 



7 (D.15) 



where 770 is rewritten as 77(A). For the space dimensions, d = A — e, the angle 
integral is given by 

where is the surface area of a unit sphere in 4 dimensions divided by (Stt)^. 
The RG equation in 4 — e dimensions reads 

-A^ = ^/(A)A(A), (D.18) 

where we introduced 

for a later convenience. 

By making coarse graining of the viscous mode with a similar procedures as 
above, we obtain a small correction for the shear viscosity: 

(5r/ = - hm [(fc2/i-i(d-l))-i V(E^^(fe,c.)y, (D.20) 

i 

where (Sxx(fc, a;))ij is the self energy for the viscous mode and given by 

Jq 

X x^'_(g)-xo'(fc-g) _ P 21) 

-iuj + Ao(9)xo ^(9) + Kk - q)Xa ^(k-q)' 
which is diagrammatically represented as Fig. ID. 41 . In the space dimension 
d = A — e, wc find the RG equation for the shear viscosity 

where the prefactor 1/24 comes from the angular integral in Eq. (|D.2ip and the 
factor {d - 1)-^ in Eq. jElQl). 

Differentiating Eq. (|D.19|) about A, we have the RG equation for /(A) 

-A^ = (6-g/(A))/(A). (D.23) 



66 



Appendix E 

Renormalization of the 
sound mode 



Next, let us make a coarse graining of the sound mode for the renormahzed 
bulk viscosity. Because a feedback from the sound mode is neglected, we must 
renormalize the mode with a method separating relevant and secondary modes 
[55] . Here, we take the method developed by Onuki [371 121] , in which RG 
equations are derived from fluctuation-dissipation relations. 

Now, we consider the equation of motion for the sound mode, (|C.22[) : 



— iui5J\\ 



-ikTr 



where the noise term 0jj satisfies the fluctuation dissipation relation: 

{e\{ki,ui)~e\{k2,uj2)) =2T,fc2^^ 



(E.l) 



(E.2) 



Since 5J\\ is a conserved density projected onto fc, we can rewrite Eq. (jE.l[) as 



iujSJn {k, uj) ^ ik ■ n(fe, ijj) ■ k, 



(E.3) 



where Hij is the stress tensor. If we take z direction as fc, Tl^z reads 



fl,,(fc,w) =-T, 



XQ\k)Cjijj{k,uj) +Co^l32hcrh{k,uj) 



+^2 hc-fo I i^{q)tp{k - q) 

qfl 



-ikiy^h^ (5J||(A;,a;) + 7rj^(fe, oj), 
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(E.4) 
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where 7f^^ (fe,aj) is the random-stress tensor coming from microscopic process 
and satisfies the relation, ik ■ 7r°(fc,w) • k ~ 9'^^{k,uj). 

We now consider how Eq. (|E.ip is affected by the coarse-graining procedure. 

~ ~ s 

In the coarse-graining procedure, the variables, ip^ , rh^ and SJ are eliminated 



from Eq. (jE.l[) . The eliminated variables do not disappear from the equation 
of motion but are implicitly contained in the noise term. In other words, we 
convert the macroscopic process in the wavenumber shell A — 5A < k < A 
into the microscopic process. In this procedure, the noise term is implicitly 
renormalized as follows 



where 



~9f{k,w)^^{k,iu) + 9f^'''°{k,iu), (E.5) 



5i[l^^"°(fc,w) =ik-n^^"°{k,uj)-k, {E.> 



-Macro / 7^ \ rri 



+ ikv\^h-^5Jl, (E.7) 



+P2 hclo / r{qW{k-q) 
Jqn 



IqQ. 

where we neglect the linear terms in Eq. (|E.7|) that is irrelevant for the following 
argument. The new term 6'||^'*'^''°(fc, a;), being due to the coarse graining, con- 
tributes the transport coefficient through the fluctuation-dissipation relation: 

= 2Tckl5v\ki,u:{){2TiY+^5{ki + fca), (E.9) 

where we have assumed that the renormalized equation of motion has the same 
form as Eq. (IE.1[) . We note that this assumption is equivalent to the requirement 
below Eq. (|D.10p . Now, we calculate the left-hand side in Eq. (IE.9|) : 

X / (^^(gi)Vi^(fci~<zi)V;S(g2)^"(fc2-g2)). (E.IO) 
Approximating the variable by the bare one, ^ ij;^^^ we find 

( 0~^-™(fci,C.i)0~^-™(fe2,L.2)) = (27r)'^+i<5(fci + fc2) 



x2fcf(T,M2"')'7o / C;^iq)Cll^iki~q), (E.ll) 

where we have used Eq. (jC.62[) and neglected a term corresponding to a dis- 
connected diagram. Then, comparing with Eq. (|E.9p . we obtain the correction 
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to the longitudinal-kinetic viscosity: 

6u\k,uj) = T^P^^hl^l [ Cf{q)Cf{k - q). (E.12) 

We are not interested in the frequency- or wavenumber-dependent bulk viscosity 
and then take the limit — >^ 0: 

Sv'' = lim 5i/\k,Lo) 

= T,P^^hhl [ (c^J'iq))'. (E.13) 
Jqn \ / 

After the integration, we find the RG equation for longitudinal kinetic viscosity: 

where we have rewritten the static parameter 70 as 7(A) to denote its cutoff 
dependence as mentioned in the text. The asymptotic behavior obtained from 
this RG equation is different from the shear viscosity's behavior, so we replace 
above RG equation as 

Although, by this method, we could more easily obtain the RG equations for 
the thermal conductivity and shear viscosity, we have taken the diagrammatic 
method for an instructive purpose. 
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